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^ . Abstract 

^vq I In this talk I first briefly explain the puzzle presented by the high muon 

content observed in air-showers which point to the ultra high-energy 7 sources 
. , like the HER-X-1. Since it had been suggested that a possible explanation of 

(y-^ I the puzzle might come from the effects of the hadronic structure of the photon, 

^D ■ I briefly explain the concept of photon structure function and comment on the 

^^ . uncertainties in the predictions of the cr^i. Then I show that while the current 

f^ I experiments at the e+e~ (TRISTAN and LEP) and ep (HERA) colliders have 

^' [ seen clear evidence for the hadronic structure of the photon, the observed muon 

excess in the air shower experiments, if confirmed by other experiments to be 
at the same high level, can not be explained in terms of the photon structure 
function. 



a^ 



Ph: 



1 What is fi puzzle? 

i^ ! One of the aims of the TeV/PeV 7 ray astronomy is to gain information about the 

c^ I origin of cosmic rays |]1[ . The extensive air-shower array experiments look for point 

sources of ultrahigh energy 7 rays p|, |^, ^, ^. These experiments can serve as di- 
rectional cosmic 7 ray telescopes. But they have to be able to distinguish between 
photon and hadron initiated air-showers. One of the criteria that is normally used is 
the expected low /i content of the 7 showers. The 7 initiated showers are supposed to 
be fi poor and the reasoning goes as follows. The muons in the air-showers come from 
TT — > /iz/ as well as from the Bethe-Heitler production of /i pairs and at still higher 
energies from the heavy quark decays. In case of 7 induced air-showers, the total 
cross-section is dominated by the Bethe-Heitler production of e~^e~ pairs. Hence one 
expects the /i content of the 7 initiated showers to be at few % level that of (roughly 
a factor 30 below) the hadron initiated showers. 
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However, so far there never has been evidence for directional showers associated 
with point sources of 7 rays which are /z poor. On the other hand, the experiments 
reported air showers associated with known 'point' 7 sources whose /x content is 
consistent with that of the hadron initiated showers or even more. Ref. reported 
observation of air showers associated with the point source Cyg X-3, ref. reported 
air-showers associated with Crab-Nebula. Air showers reported in ref. |Q were as- 
sociated with the point source HER X-1 which were /i rich, i.e., their n content was 
consistent with that of a hadronic shower. This observation as well as any failure 
to see a fi poor shower associated with 'point' sources, signalled existence of what is 
termed as the /z puzzle, i.e., the cosmic ray air shower experiments 'see' more muons 
than they 'ought' to. 

Let us now turn more specifically to the theoretical predictions H, |^, |, § of the 
H content of the 7 induced showers. The actual predictions involve detailed Monte 
Carlo simulation of the shower developement, but the essential features of how one 
arrives at these predictions can be summarised as follows (see, e.g., ref. 0). One 
starts with measured photoproduction cross-sections in the laboratory. In the days 
before the results from the ep collider HERA at DESY, the maximum centre of mass 
energy {\^)^p for which data were available was ~ 20 GeV. At these energies the 
atyp is ~ 100 fib. Now to estimate the /i content of the airshowers one has to 
extrapolate these cross-sections to TeV/PeV energies using the experimentally well 
known logarithmic rise of total cross-sections. Using this one then calculates the 
production of hadrons and hence production of muons, telling us how many /i's ought 
to be there in the 7 induced showers. The observed /i excess [Q would require a ai^t, 
~ 0(100) mb for PeV energy photons. This would imply a faster than the assumed 
logarithmic growth of a^p in going from the laboratory measurements at GeV energies 
to the PeV energies involved in the cosmic ray experiments. This can happen only 
if there exist some new threshold in the photonuclear (7 air ) cross-sections which 
gives it a steeper energy dependence. Hadronic structure of the high energy 7 (in 
particular the gluon content) can indeed cause a sharp increase in the cross-sections 
with increasing 7 energy [|l^ . It should be noted here that the increasing importance 
of the hadronic structure of the 7 with rising energy was not introduced here to 
'explain' the // 'puzzle', but is a prediction of perturbative QCD. What is not clear, 
and hence has been a topic of debate, is how much does the hadronic structure of 
photon contribute to the rise in a^t, with energy. 

2 Hadronic content of photon 

The terminology of the 'structure ' of a photon is essentially a short hand way of 
describing how a high energy photon interacts with other particles: hadrons and 
photons. The idea that photons behave like hadrons when interacting with other 
hadrons dates back to the early days of strong interaction physics and is known to 
us under the name of the Vector Meson Dominance (VMD) picture. This essentially 



means that at low 4-niomentuni transfer, the interaction of a photon with hadrons is 
dominated by the exchange of vector mesons which have the same quantum numbers 
as the photon. While this picture works reasonably well for 'soft' processes (i.e., 
reactions characterized by small 4-momentum transfer), it is not at all clear that 
it should describe the whole story of interactions of photons with hadrons at high 
energies as well. In the VMD picture one then expects 

„tot „, „ ^tot 

a^p ocaavp, 

where a is the fine structure constant. However, since the photon 'behaves' like a 
hadron while interacting with other hadrons it must be possible to get information 
about the photon structure just like the other hadrons, e.g., the proton. This in- 
formation is obtained by studying the deep inelastic scattering (DIS) of high energy 
leptons of energy E off proton targets. 
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The double differential cross-section for the process is a function of two independent 
variables y = u/E where u is the energy carried by the probing photon in the labora- 
tory frame, and x = Q^ /{2Mv) where M is the proton mass and —Q^ is the invariant 
mass of the virtual photon in fig. |l|a. In the quark-parton-model (QPM) it is given 



by, 



where 



Figure 1: Deep Inelastic Scattering for the proton and photon. 
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are the two electromagnetic structure functions of the proton (in the QPM Fl{x) is 
identically zero but not so in QCD) and fq/p{x) the probability for quark q to carry 



a moinentuin fraction x of the proton and Cg denotes the electromagnetic charge of 
quark q in units of the proton charge. 

To measure the structure function of a photon such an experimental situation is 
provided at e~^e~ colliders in 7*7 reactions as shown in fig. |1] (b). Here the virtual 
photon with invariant mass square —Q^ probes the structure of the real photon. If the 
VMD picture were the whole story then one would expect that such an experiment 
will find 

Then with increasing Q^, the structure function F2 will behave just like a hadronic 
proton structure function and shrink to lower values of x as predicted by QCD |TT 



However, there is a very important difference in case of photons, z.e., photons possess 
pointlike couplings to quarks. This has interesting implications for 7*7 interactions 
as first noted in the framework of the QPM by Walsh [|12[ . It essentially means that 
7*7 scattering in fig. |l] contains two contributions as shown in fig. 0. The contribution 



Figure 2: Two contributions to F2 . 

of fig. 1^ (a) can be estimated by eq.(|]), whereas that of fig. || (b) was calculated in 
the QPM [|1^]. This is done by considering the cross-section for the reaction 

7 + 7* -> g + g. 

Due to t and u channel poles this can be calculated only when one considers quarks 
with finite masses. The result can be recast in a form equivalent to eq. (0): 

= 7-^ X — 

dxdy Q'^ 71 

( 1^2 

Y: eU (1 + (1 - y?) X [x{x' + (1 - x)') X In — 
T I S 

+8x'^{l-x)-x]-y'^[4x'^{l-x)]\, (4) 

where W'^ = Q^{1 — x)/x. On comparing eqs.(0) and (Q) we see that the factors 
in square brackets in the above equation have the natural interpretation as photon 



structure functions F2 and Fl and one has 
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Two points are worth noting: the function ^^'P°™*' '^(x,Q'^) can be completely cal- 
culated in QED and secondly this contribution to F2 increases logarithmically with 
Q^. So in this simple 'VMD + QPM' picture, F2 consists of two parts, i^^'P°™ 
and F2' , with distinctly different Q^ behaviour and with the distinction that for 
one part both the x and the Q^ dependence can be calculated completely from first 
principles. 

This QPM prediction received further support when it was shown by Witten 

13| that at large Q^ and at large x, both the x and Q^ dependence of the quark and 

gluon densities in the photon can be predicted completely even after QCD radiation is 

included. An alternative way of understanding this result is to consider the evolution 



equations |T^ for the quark and gluon densities inside the photon. These contain an 
inhomogeneous term on the r.h.s proportional to a, which describes 'J —^ qq splitting, 
i.e. the pointlike coupling of photons to quarks. In the 'asymptotic' limit of large Q^ 
and large x, the fq-/^{x, Q^) have the form 

f-J-^{x,Q')^axln(-^)F,{x) 

^ —F,{x), (6) 

O-s 

where Aqcd is the usual QCD scale parameter, as(Q^) is given in terms of the running 
strong coupling constant by ^°^^ ' and the x dependence of the -Fj(x) is completely 



calculable. Note here the factor In ( -j^ — j on the r.h.s. Measurements |T^ of the 

photon structure function F^ in 7*7 processes did indeed confirm the basic QCD 
predictions of the linear rise of F2 with In (Q^) at large x. This discussion thus 
means that just like one can 'pull' quarks and gluons out of a proton one can look 
upon the photon as a source of partons and that the parton content of the photon 
rises with its energy. Physically this means that the photon splits in a qq pair and 
these radiate further gluons and thus fill up a volume around photon with partons. 

The asymptotic solutions discussed above, though very useful to understand the 
rise of the photon structure function with Q^ , are valid only at large x and large Q^. 
At small values of x these solutions diverge, indicating thereby that 'hadronic' part 
of F2 can not be neglected at small x. Hence in practice, this separation of F2 in 
two parts as in fig. |^ is not very meaningful, especially when one wants to use this 
parton language to predict the high energy photon interactions. Although the debate 



on the subject is not yet closed |TB[ , it is now generally accepted that it is better 
to forego the absolute predictions of F2 of the asymptotic part, that are possible in 
perturbative QCD (pQCD) and use only the prediction of the Q^ evolution of the 
photon structure function in analogy to the case of the proton structure function. At 
present there exist eight different parametrisations of the photon structure function 
|I^ |T^. The DIS measurements described above measure only the quark-parton den- 
sities fq^/^{x, Q2) (for X > 0.05 and Q^ < 100 - 200 GeV^) directly and fg/^{x, Q^) is 
only inferred indirectly. As a result there is considerable uncertainty in the knowledge 
of fg/.y{x, Q"^). The different parametrisations differ quite a lot from each other in the 
gluon content. It should also be mentioned here, that these differences reflect the 
differences in different physical assumptions in getting /^/^(a;, Q^) from the data on 
F2 . So independent information on fg/^^XjQ"^) is welcome. 

3 Calculation of jet production in 77, ^p collisions 

In this section let us discuss how one can compute the high energy 7 cross-sections 
using the parton language and what are the crucial factors affecting these predictions. 
In the parton language, interaction of a hadron with others can be described, at high 
energies and for processes involving final state particles at large angles to the original 
beam direction (large transverse momentum px), in terms of the scattering of the 
pointlike constituents inside the two hadrons against each other. In this picture, the 
differntial cross-section for the production of a pair of two large pt jets in the collision 
of a 7 with a proton(say) will be given by 



da 



{.IV -^ jets + X) = Y^ j dx^ fpif-ri^i) j dxp / 
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^P^ Pl,P2,P3,P4' 

X^{P,+P2^P,+P,), (7) 

dpT 

where the sum is over all the different intital (final) state partons Pi,P2 {P3,Pa) 
and j^ is the subprocess cross-section that can be computed in pQCD. As we have 
already seen in the discussions on F2 one can get effectively 'real' high energy 7 
beams in the laboratory by using e-beams and then making sure that the final state 
e is scattered in the forward direction at a very small angle. In this situation the 
large pr jet production in high energy e — p collisions can be computed as 



— (ep -^ jets + X) = Yl J dz f^/e(,z) J dx^ fp^^ 



da , 
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X j dXp fp,/p{Xp) -^ (Pi + P2 -^ ^3 + P4) ■ 



The 7 induced processes are of two types : 



Figure 3: 'Direct' and 'resolved' contributions to jet production in ep collision. 

1) The 'direct' processes, an example of which is depicted in fig. ^ (a), where the 7 
couples directly to the partons in the photon. In this case all the energy of the 
photons goes into the subprocess and and the hadronic content of the 7 plays 
no role and in eq.(|^), 

2) The 'resolved' processes where the partons in the proton interact with the par- 
tons in the photon and hence only a partial fraction of the 7 energy is available 
for the subprocess. One of the possible contribution is shown in fig. (b). 

It should be noted here that although the 'resolved' processes have an extra factor of 

as, due to the factor of — in the parton desities of the photon (recall eq. (Bf)), both 

as 
the 'direct' and the 'resolved' processes are formally of the same order in the coupling 

constants. With rising 7 energies, increasingly more energy becomes available for 
the subprocess, at a fixed pt or mass of the final state. Hence the importance of 
the 'resolved' processes increases with the 7 energies. The gluon content of photon 
/g/^(x, Q^) is peaked at small values of x and hence the importance of the contri- 
butions involving the gluon in the photon in the initial state in eq. (H), increases 
with the increasing 7 energy at a fixed pt or with decreasing p^ at a fixed 7 energy. 
The 'resolved' events will also have additional 'spectator' jets in the direction of the 
incident 7 {i.e. the direction of the incident e). 

Such high photon energies are available at the HERA collider at DESY (Hamburg) 
in the collision of a 30 GeV e beam with a p beam of 800 GeV. This corresponds to 
cm. energies for the jp system < 300 GeV (which in turn means E^ < 40 TeV in 
the frame where the proton is at rest). Indeed a calculation [^ showed that at the 
HERA collider the photo-production of jets in the process 

7 + P ^ jets + X, 

is dominated by the 'resolved' processes upto px — 40 GeV. This dominance is sen- 
sitive to the gluon content of the photon and fortunately not very sensitive to the 
choice of gluon parametrisation of the proton, partially because of the much better 



Figure 4: Ratio of resolved and direct contributions for da{ep -^ iets) /dp^ as a 



function of p-r |18]. For details see ref. [18 



knowledge of fg/p in the relevant Xp region. An example is shown in fig. ^ for two dif- 
ferent paranietrisations of fgi^{x,Q'^) that were then availble. Since then the HERA 
experiments HI and ZEUS |T^, ^ have studied the photo-production of the jets, and 
confirmed the existence of the 'resolved' contribution at the expected level, verified 
various expected qualitative features of the resolved contributions such as existence of 
the spectator jets in the backward direction, different angular distribution expected 
for the jets produced from the hard scattering of the partons in the photon etc. An 
example of the same is shown in fig. |^. As a matter of fact the 'resolved' contributions 



Figure 5: The total ep cross-section measured [^ for transverse energies larger than 
E^. The curve is the HERWIG prediction, using the DG parametrization with pT,min 
= 1.5 GeV. 



to the jet-production have been isolated and used to extract fg/^y^x, Q"^) |^ 



Similar studies of the jet-production in 77 collisions [22| at the e^e collider TRIS 



TAN and LEP, have confirmed the existence of the 'resolved' contributions vM- These 



studies have already ruled out some of the very hard parametrisations of fg/-y{x, Q^) 
p2 , 24]. Thus these observations have provided a confirmation (in addition to the 



DIS measurements ) of the ideas about F2 and these experiments will continue to 
add to our knowledge of the fg/-y{x, Q^), fq^^^x, Q"^). 

4 F2 and QCD prediction of atSL 

What is more relevant for the issue of ft puzzle is the total photopro duct ion cross- 
section cr^p. The calculation of the 'hard' processes such as the jet-production cross- 
sections is well defined in pQCD but valid upto ~ 1-2 GeV. The total inclusive jp 
cross-section for production of jets with pt > p™"^ given by, 

<., = I ^{IP- Jl +J2+X) dpT, (9) 

rises very strongly with decreasing pT,min at a fixed 7 energy or with 7 energy at a 
fixed PT,min- The differential cross-section ^(7P — > J1+J2+X ) receives contribution 
from the 'direct' as well as the 'resolved' processes as discussed before and it also rises 
strongly with decreasing pt- At high jp cm. energies the resolved photon processes 
then cause copious production of the 'minjets' and a very rapid rise of the inclusive 



7P cross-section, as was first pointed out in ref. [jTO|- An example is shown in fig. 



Figure 6: Predictions [|1^ of the increase of the inclusive (mini)jet cross-section in 
7P collisions with y/s, for pT,min = 2 GeV and various parametrizations for f^/^. 

for pT,min = 2 GeV and various parametrizations of f^/^. Of course, the total cross- 
section cannot grow indefinitely at the rate shown in fig. |^; some mechanism will 



have to unitarize it. This problem is well known for hadronic {pp or pp) collisions. In 
this case unitarization is usually achieved by eikonalization. The crucial observation 
here is that LO QCD predictions for cross-sections, like those shown in fig. P, refer 
to inclusive jet cross-sections; in other words, they differ from the jet production 
contribution to the total cross-section by a factor of the average jet pair multiplicity 
(njet). Formally one writes (following ref. |2B] and modifying for the 7p case as 
pointed out in 



<-i = J d% Phad {l-exp [- {a];f{s) + xt) Mb)/P^ } , (10) 

Here b is the two-component impact parameter, A{b) describes the transverse 
overlap of partons in the nucleon and the photon, a^^p^"^ is the perturbative QCD 
prediction for the minijet cross-section (obtained by integrating da/dpr in the region 
Pt > PT,min), and x^pp^ is the non-perturbative (soft) contribution to the eikonal, 
which is fitted from low-energy data and Phad is a parameter describing the proba- 
bility that the photon goes into a hadronic state; clearly Phad ~ 0{a). Thus we see 
that unlike the predictions of the jet cross-sections the predictions of a^p depend 
not only on the fg/^{x,Q'^) and fg./^{x,Q'^) but also on pT,min, Phad, A{b) and the 
nonperturbative contribution to the eikonal. There is considerable theoretical uncer- 
tainty in the ansatze uses for cr'^"^^ and the choice of Phad, as well as PT,min and A{b) 
Thus the predictions of the total 7p cross-sections will depend on all 



these parameters apart from the information about fg/^{x,Q ) and fq,/-^{x,Q ). 

Thus to recapitulate the predictions for at^t, in the frameowrk of QCD depend on 
the following factors: 

1. The parton densities in the photon, particularly fg/^^x^Q"^). The range of x 
values that are relevant shifts to smaller values from x < 0.01 to x < 0.001 or 
smaller, as one goes from HERA energies to the case of the PeV energy photons 
relevant for the fi puzzle. 

2. Value oi pT,min- 

3. Modelling of the soft cross-section. 

4. The probability Phad for the proton to go to a hadron state. 

5. The transverse overlap A{b) of the partons in the nucleon and the photon. 



Fortunately the experiments at HERA [30, 31 1 also measured a!^L in addition to 



HP 



electro-(photo-)production of 'hard' jets. This information from HERA when coupled 
with information about the multi-jet production in 77 collisions p2|, |2^, can be used 
to reduce the uncertainty in the knowledge of /^/^(a:, Q^) and PT,min- The 77 data 
are consistent with pT,min — 2.0 — 2.5 GeV depending upon the choice of parametrisa- 
tion for fg/^i^x^Q"^) and fqi/'y{x,Q'^). But the very broad gluon distribution of LACS 



parametrisation is ruled out already by these data. The observation of the resolved 
contribution to large pr jets at HERA [|1^, |2y, ^ allows an extraction of fg/-y{x, Q"^) 



which is essentially consistent with almost all the modern parametrisations (that rise 
somewhat steeply at small x) other than LACS and a similar value of PT,min- It 
should also be mentioned here that the present 'extractions' of /g/-y(a;, Q^) has large 
experimental as well as (Monte Carlo related) theoretical errors. Both these errors 
are likely to go down in further analysis/data from HERA, allowing perhaps a better 
discrimination. 

As seen above, calculation of a^L involves, in addition to pT,min and fg/^{x,Q'^), 
fg./^{x,Q'^), Phad, A{b) as well as model for a^°^^. Here some model-builders take 
recourse to the low-energy measurements of a!^p to determine the model parameters. 
Here almost the whole cross-section is dominated by the soft process and the only 
relevant 'hard' contribution is the 'direct ' process. As a result of these the low energy 
predictions are independent of /^y^ and can determine some of the nonperturbative 
parameters, especially x^°^^ quite well. The ZEUS measurement |^0| along with the 
(pre-HERA) predictions using the eikonalised calculation mentioned above for differ- 
ent choices of the parameters mentioned are shown in fig. |^. The two solid curves show 



Figure 7: Comparison of various predictions of total ■yp cross-sections with low-energy 
data and the recent ZEUS measurement |^ . 



fits to low-energy data based on Pomeron phenomenology which does not involve any 
'hard' contribution. The two dot-dashed curves show minijet predictions [^ using 



the DC parametrization with pT,min = 1-4 (upper) and 2.0 (lower curve) GeV, while 
the dotted and dashed curves have been obtained from the LACl parametrization 
using the same values of PT.min- The LAC parametrization seems disfavoured but 
in view of the theoretical and experimental uncertainties it might be premature to 
exclude it altogether. The DC minijet prediction with pr,mm = 2 GeV is certainly 



in agreement with the data. At the HERA energy the effect of eikonahsation on the 
predictions for a^p is not very large, but it is certainly substantial when the PeV 
energy 7 are involved. A further clarification between different models (eikonalised 
minijet or Pomeron based) will be possible if one can measure energy-dependence of 
a^p. At present, one can say definitely: 

1. The resolved contribution to photoproduction of jets does rise with increasing 
7 energies and has been seen. 

2. These processes also cause the a^L to rise with energy. However, with pT,min — 2 
GeV and almost all the current parametrisations of /^/^(a;, Q"^) and fg./^{x, Q^) 
this rise is much less steep than the most 'optimistic' pre-HERA predictions of 

tot 

u^p. 

3. Further measurements of jet cross-sections, details of event shapes, e.g, mul- 



tiplicity and various correlations ||3^ should provide further help in clarifying 
this situation. 



tot 



Armed with this information now one can turn towards the predictions of a^p 
for the TeV/PeV energy photons from the point 7 sources. Fig. ^ shows predictions 



Figure 8: cr^^^l^ and crj^gj cross-section as a function of 7 energy in an eikonalised 
model. lESI . 



of a particular model ||2^ for aj^er. This model of •yp interactions fits the PT,min 
from IT N scattering data and the quark-gluon densities of the 'hadronic states' into 



which the photon converts itself are related to parton-distributions for a pion. The 
model predictions, with the inclusion of soft scattering background, agree with data 
at low-energy (8 GeV < \^ < 20 GeV ) quite well and are reasonable (though 
somewhat high) for HERA energies. As can be seen from this figure o"inei does rise 
with increasing E^ (or equivalently y/s^p) and is about ~ 2 — 3 mb for the PeV 
energy photons. This falls way short of the a ~ 100 mb required to explain the muon 
excess observed in the CYGNUS experiment. One should keep in mind that since 
(rijet) does rise sharply with energy (recall fig. |^) it is possible that fluctuations could 
sometimes produce /i rich showers. But if the future air shower experiments continue 
to see showers associated with point sources, with /i content as high as that seen by 
the earlier experiments , the explanation must lie somewhere else and not with the 
rising 7p cross-sections 

So in conclusion we can say that 

1. There exists a hint of a /i excess seen by the extensive air shower experiments 
in the air showers associated with 'point' sources. This would mean that the 7 
induced showers look more like hadronic showers. 

2. Since the photon does behave like a hadron at high energies and the 'resolved' 
contribution of the 'partons in the photon' increases with 7 energy, feature (1) 
could have been, in principle, due to the hadronic content of 7. 

3. The current collider experiments at e~^e~ (effective 77) colliders TRISTAN and 
LEP and ep (effective 7p) collider HERA have seen these contribiutions coming 
from the hadronic structure of the photon and have provided some information 
about fg/^{x,Q^). 

4. The prediction of a^L has considerable theoretical uncertainties, but using the 
current experimental information from the e+e~, ep colliders on jet production 
and on a^p from HERA, one can conclude that the hadronic structure of 7 can 
cause an increase in a!^p by a factor of ~ 2-3 for the PeV energy photons but 
not much more. 

5. Hence the photon strucure function effects can give rise to n rich showers only as 
a fluctuation but if future experiments continue to report existence of extremely 
H rich showers , the explanation can not be provided in terms of the photon 
structure. 
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